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Abstract. Quantum phase transitions (QPTs) and the ground-state phase diagram 
of the spin-1/2 Heisenberg-Ising alternating chain (HIAC) with uniform Dzyaloshinskii- 
Moriya (DM) interaction are investigated by a matrix-product-state (MPS) method. 
By calculating the odd- and even-string order parameters, we recognize two kinds 
of Haldane phases, i.e., the odd- and even-Haldane phases. Furthermore, doubly 
degenerate entanglement spectra on odd and even bonds are observed in odd- and even- 
Haldane phases, respectively. A rich phase diagram including four different phases, 
i.e., an antiferromagnetic (AF), AF stripe, odd- and even-Haldane phases, is obtained. 
These phases are found to be separated by continuous QPTs: the topological QPT 
between the odd- and even-Haldane phases is verified to be continuous and corresponds 
to conformal field theory with central charge c=l; while the rest phase transitions in 
the phase diagram are found to be c=l/2. We also revisit, with our MPS method, the 
exactly solvable case of HIAC model with DM interactions only on odd bonds, and 
find that the even-Haldane phase disappears, but the other three phases, i.e., the AF, 
AF stripe, and odd-Haldane phases, still remain in the phase diagram. We exhibit the 
evolution of the even-Haldane phase by tuning the DM interactions on the even bonds 
gradually. 
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It is well known that the gronnd-state properties are intrinsically different 
between the nniform spin-1/2 and spin-1 Heisenberg antiferromagnetic chains (HAFCs). 
Specihcally, the ground state of the former (spin-1/2) model is critical and has gapless 
low-energy excitations, he., spinons with fractional spin pQ; while for the latter spin-1 
HAFC, there exists a hnite energy gap, he., the Haldane gap, above the ground state, 
and the elementary excitations are magnons with nonzero mass[2]. This intriguing 
feature of the spin-1 HAFC was hrst conjectured by Haldane |3], and then theoretically 
interpretated in terms of the valence bond solid state picture [1]. In addition, a 
nonlocal string order parameter Og has been introduced to characterize the “dilute” 
antiferromagnetic (AF) order in the Haldane phase [S]. Further study shows that the 
topological long-range order in the SO{2n + 1) symmetric matrix product state (MPS) 
can be fully identihed and characterized by a set of nonlocal string order parameters [0] . 
The SO{2n+l) symmetric MPS contains diluted AF orders in n different channels and a 
hidden (Z 2 XZ 2 )” symmetry breaking, which opens an excitation gap. The existence of a 
nonlocal order parameter and a dilute AF order actually are physically observables: the 
experimental observation of spin-1/2 edge states in spin-1 chains evidences the valence 
bond picture in the Haldane phase miH], and it is worthy noticing that people recently 
succeed to observe directly the nonzero string order in cold atom systems [Hj. 

According to Landau’s phase transition theory mm, different phases and their 
phase transitions can be characterized by local order parameters reflecting different 
symmetry breakings. However, the Haldane phase, which has no explicit symmetry 
breaking, is separated from a ’’trivial” phase, e.g., large-H phase {D represents single¬ 
ion anisotropy), by a phase transition [121 ESI E], as long as symmetries [SO(3), bond- 
centered inversion, etc] are intact in the Hamiltonian space (path). So to speak, the 
Haldane phase has a symmetry-protected topological (SPT) order [15], and the string 
order parameter Og can thus be regarded as a quantihcation of the topological order of 
Haldane phase. Moreover, a nonlocal unitary transformation was constructed to uncover 
explicitly the hidden Z 2 x Z 2 symmetry breaking in the Haldane phase [Ml El EHl [I9|. 

On the other hand, SPT phases have also been reported in spin-1/2 multi-period 
chains. A spin-1/2 ferromagnetic-antiferromagnetic alternating chain has been shown to 
have nonzero string order, which changes continuously through the Haldane and dimer 
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phases and show a crossover behavior between them Bni. In addition, Kohmoto and 
Tasaki indicated that the spin-1/2 ferromagnetic-antiferromagnetic alternating chain 
fnlly breaks the hidden Z 2 XZ 2 symmetry not only in the strongly coupled (Haldane) 
phase but also in the decoupled (dimer) phase, supporting that the spin-1/2 dimer state 
belongs to a Haldane-like phase. Recently, Wang et al [21] associates the phase transition 
between an odd- and even-dimer states as topological quantum phase transition (TQPT) 
between two kinds of Haldane phases (even- and odd-Haldane), which have different 
string order parameters. 

Among others, the Heisenberg-Ising alternating chain (HIAC), originally proposed 
by Lieb et al 122] and re-examined subsequently by Yao et al [23], constitutes an 
interesting spin-1/2 chain model. The Hamiltonian of the HIAC is given by 

Af/2 

H = Y,IJhSx-i ■ S 2 . + (1) 

i 

where S denotes a spin-1/2 operator is the 2 ;-component), and N (an even number) is 
the total number of spins. Jh and Jj {Jh, Ji > 0) are the Heisenberg and Ising couplings 
on the odd and even bonds, respectively. This model can be solved exactly, and a phase 
transition from a quantum paramagnetic phase to the AF phase was found to occur at 
Ji/Jh = 2. When the Dzyaloshinskii-Moriya (DM) interactions are switched on odd 
bonds of the HIAC [23], the model is still exactly solvable, it turns out that the ground- 
state energy exhibits an interesting non-analytic behavior accompanied by a gapless 
excitation spectrum along the line Ji=2\jD'^ + Jjj, where D denotes the strength of 
DM interaction [see Eq. ([^ below]. This critical line was found to separate a quantum 
paramagnetic phase (J/ < 2^/0“^ + Jfj) and an AF one. However, in both cases, the 
property of the SPT order has not been discussed in the paramagnetic phase. 

In this paper, we study an extended HIAC model, which has an uniform DM 
interactions on both even- and odd-bonds [Eq. ([^ below], which is thus no longer exactly 
soluble. We here adopt a MPS based numerical method, i.e., the infinite time-evolving 
block decimation (iTEBD) algorithm[25l [26], to determine the ground state with high 
accuracy. We obtain a rich ground-state phase diagram which includes an AF phase, an 
antiferromagnetic stripe phase (AFSP), and two Haldane phases. In particular, the two 
Haldane phases, namely, the odd- and even-Haldane phases, can be characterized (and 
distinguished) by two different nonlocal string order parameters and doubly degenerate 
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entanglement spectra on the odd and even bonds, respectively. Moreover, a continuous 
TQPT between the odd- and even-Haldane phases has been found. 

The rest of this paper are arranged as follows. In section]^ the model Hamiltonian 
and the exact phase diagram is present. In section the numerical results are shown 
and analyzed in details. Lastly, we devote section]^ to a discussion and summary. 

2. Model Hamiltonian and the phase diagram 

The spin-1/2 HIAC with an uniform DM interaction is described by 

N/2 

H = . S2. + - Sjt-iSS) 

+ + £>(51.51,+, - S».Sf.+i)]. (2) 

Here D is called the DM vector, with D=Dez adopted in this paper. The DM 
interaction is an antisymmetric spin-spin coupling, due to the spin-orbit coupling effect. 
Recently such DM-like spin-orbit interaction in the single-crystal yttrium iron garnet 
was experimentally verified and found amenable to manipulations [27] . It is worth noting 
the DM interaction can be eliminated by a canonical spin rotation pTl 128] , and then 
Eq. ([^ can be mapped into a period-two XXZ model consequently. In contrast to Eq. 
Q, the transformed XXZ model is no longer exactly soluble. 

We employ the iTEBD method to investigate the ground-state properties of model 
Eq. (|^. The iTEBD algorithm is a numerical approach based on the MBS ansatz, 
which accurately describes quantum many-body wavefunctions in one dimension, given 
sufficiently large bond dimensions y. During iTEBD optimizations, one successively 
applies imaginary-time evolution gates exp {—rh) on a trial MBS state, until the latter 
converges to the variational ground state, h is the local Hamiltonian (here only two-site 
coupling term on an even or odd bond), and r is the Trotter step length. In practice, 
we start from r = 10“^ and gradually reduce it to smaller values (eventually down to 
r = 10“®), with a total evolution steps of number O(10^~^), given the variational MBS 
wavefunction is well converged in the course of imaginary-time evolutions. Regarding 
the bond dimension y; of MBS used in practice, we found for this specific model, y ~ 30 
can already provide rather accurate and converged results (for various gapped phases of 
the model). Remarkably, distinct from other finite-size algorithms like quantum Monte 
Carlo or finite-size density matrix renormalization group method, iTEBD exploits the 
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translation invariance of the one-dimenional chain and make the thermodynamic limit 
directly accessible. 

From the numerical results, we conclude a rich ground-state phase diagram of 
the spin-1/2 HIAC model with DM interactions, which is shown in Fig. One can 
find four different phases: an AFSP, AF, and two distinct Haldane (odd- and even- 
Haldane) phases, all separated by critical lines. Two magnetic ordered phases can be 
identified by local order parameters (stripe and Neel orders, respectively), while the non¬ 
magnetic odd- and even-Haldane phases are distinguished by two nonzero string order 
parameters. In Fig. besides the QPTs between different magnetic ordered phases and 
those between ordered and Haldane phases, we also observe a TQPT between the odd- 
and even-Haldane phases, which can not be explained via symmetry breaking (Landau 
paradigm). This TQPT can be well captured by the singular behavior of the bipartite 
entanglement, as well as even- and odd- string order parameters. 

In addition, we are also interested in the critical properties on the phase boundaries. 
We select several representative points (PI to P5, see Fig. and calculate their central 
charges c by fitting the block entanglement entropies. We find c = 1/2 at the magnetic 
order-disorder transition points (PI, P2, and P3); while c = 1 at the TQPT points (P4 
and P5) between two kinds of Haldane phases . 


3. Numerical results and discussions 


In the following, we study the phase transitions along two selected paths, by fixing 
11=0.7 or 2.0 and tuning Ising couplings Jj. In subsection |3.1[ we study the line 11=0.7 
crossing four different phases (see Fig. [^, i.e., the AFSP, odd-Haldane, AF, and even- 


Haldane phases, and focus on the magnetic order-disorder QPTs. In subsection 3.2 


we focus on the path along D=2.0 line, along which the odd- and even-Haldane phases 
touch each other and a QPT was found between them. Since both the even- and odd- 
Haldane phases have the same symmetry properties (no symmetry breaking), this exotic 


phase transition between two phases can be regarded as a TQPT. in subsection |3.3[ the 
block entanglement entropies and the fitted central charges of several critical points are 
evaluated. In our calculations, the Heisenberg coupling Jh = I [see Eq. (|^] is set as an 
energy scale, and the retained bond dimension of MPS is set as y = 30. In order to check 
the reliability of our results, we reinvestigate the Heisenbeg-Ising chain with odd-bond 
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Figure 1. (Color online) Magnetic phase diagram of the spin-1/2 HIAC with uniform 
DM interactions [Eq. ([^]. It includes four different phases: an AFSP, an AF phase, 
and two kinds of Haldane (odd- and even-Haldane) phases. The filled square ■ denotes 
a tricritical point, and the vicinity of it is subject to substantial uncertainty. Five 
representative points (PI, P2, • • •, P5) on the the phase boundaries are selected, and 
the corresponding central charges will be discussed below. 


DM interactions in subsection 3.4[ and find that the disordered phase corresponds to 
the odd-Haldane phase. Next we study the effect of moderate strengths of even-bond 
DM interactions, and the evolution of the even-Haldane phase are shown pictorially in 
subsection 13.51 


3.1. The QPTs along the D=0.7 line 

In this subsection, we investigate the QPTs along the line D=0.7. The local 
magnetizations on site i, = (S^) {a={x, z}) are evaluated. We find that 

the transverse magnetization Mf vanishes completely along the line. However, local 
magnetizations M/ are nonzero in two regions: In the region of Ji < —2.54, the spin 

configuration is----|-h---or-- - —|-H-( + and - denote spin up and down 

in terms of M?, respectively), constitutes a period-four stripe-like AF ordered phase 
(dubbed as AFSP); another symmetry broken phase is between 2.0 < Jj < 2.6, where 
we observe an conventional AF Neel order, i.e., •••-I -1— ••• or ••• — I - [-■■■. 
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Figure 2. (Color online) The stripe order parameter + 

M|j_|_i)| and the Neel order parameter = J2f=ii~^YMi\ versus Ising coupling 

J/ are shown (along the D—0.7 line). 

It is convenient to define the stripe order parameter M^tripe = + 

and the Neel order parameter l)*Mf | to distingnish these 

two magnetic ordered phases. and are evaluated for various J/ and 

shown in Fig. As expected, M^^j-ipe and can be used to distinguish the 

AFSP and AF phase; however, they are not able to tell the differences between the 
rest two nonmagnetic regions (—2.54 < Jj < 2 and J/ > 2.6), where local magnetic 
moments vanish and no symmetry breaking is observed. Interestingly, by introducing 
some nonlocal order parameters, we see that these two regions (—2.54 < Jj < 2 and 
J/ > 2.6) have different types of non-local string order, suggesting the existence of two 
different SPT orders. We will discuss these two distinct phases, namely the odd- and 
even-Haldane phases, and the phase transition between them later in subsection |3.2[ 
The entanglement entropy has been shown to be an efficient tool to detect QPTs 
[2911301131]. For continuous QPTs, owing to the quantum criticality at the transition 
point, the bipartite von Neumann entanglement entropy 

5b =-Tr(plog2p), (3) 

diverges at the critical point. In the MPS framework, 5 b can be evaluated by 
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Figure 3. (Color online) Bipartite entanglement with _D=0.7 on odd bond (5'2i-i,2i) 
and even bond (S' 2 i, 2 i+i)- 

= - Ef=i , where A^’s are (normalized) singular values obtained through 

bond singular value decomposition, i.e., square root of the eigenvalues of the half- 
inhnite-chain reduced density matrix. In Fig. we show the results of the bipartite 
entanglement entropies 821 - 1 , 2 % on the odd bond and 821 , 21+1 on the even bond. 8 b 
shows three sharp peaks, indicating three continuous QPT points. Besides, we have 
also calculated the ground-state energy per site and and its derivatives (versus J/), 
and show them in Fig. Although the ground-state energy per site e* and its hrst- 
order derivative {dei/dJj) behave continuously [Fig. (a) and (b), respectively], the 
second-order derivative {d?ei/dJ‘j) exhibits divergent peaks at three critical points [Fig. 

(b)]. The singularities in d?ei/dJ‘f conhrm the conclusion that these three magnetic 
order-disorder QPTs are of second order [3^ . 

3.2. The D=2 line and the topological phase transition 

Next, we explore the line D=2 in Fig. along which the bipartite entanglement 
entropy is calculated (see Fig. [^. We will focus on the QPT occurred at Jj ~ 1.43 
(P4 point of Fig. [^, between two magnetic disordered phases, i.e., the odd- and even- 
Haldane phases. The sharp peak of the entanglement entropy curve suggests that this 
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Figure 4. (Color online) (a) Ground-state energy per site with 13=0.7; (b) The 
first- and second-order derivative curves of e^. 


QPT at P4 is also of second order. 

We firstly show the results of the nonlocal string order parameters, which is 
introduced by den Nijs and Rommels [5], and are used to characterize the dilute AF 
order in the Haldane phase of the spin-1 Heisenberg chain. The presence of a nonzero 
string order parameter can be explained via a hidden Z 2 x Z 2 symmetry breaking [BE], 
and be related to a symmetry-protected topological order. For the present model, in 
order to distinguish the odd- and even-Haldane phases, we need to dehne two different 
string orders (odd- and even-string orders), which are dehned as |32] 


O a,odd/T\ /n n n \ 

s (A) — \^2i-l’^2i' ' '1 

O a.even(T\ / \ 

s (A) — \<72i(^2i+l ■ ■ ■ ^2j^2j+l)^ 


( 4 ) 

( 5 ) 


in which L = \j —i\, and (t“’s {a=x^ y and z) denote spin-1/2 pauli matrices. The 0"’°'^'^ 
is calculated from an odd site (2i — 1) to an even site (2j), but the starts from an 

even site (2i) and ends at an odd site {2j -|- 1), i.e., the lattice distance is L=2{j — i -|- 1) 
in both cases. 

In Fig.|^ we show the results of both the odd-string and even-string 

order parameters at two representative points (Jjy=l, J/=1.2, and D=2) and (Jiy=l, 
J/=1.6, and D=2), which locate at distinct nonmagnetic phases (Fig. [^. From Fig. 
we hnd that 0|’°‘^'^(L) and show different behaviors at these two points: The 








Bipartite entanglement S 
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Figure 5. (Color online) Bipartite entanglement on odd bond ( 821 - 1 , 21 ) and even 
bond (S 2 i, 2 i+i) with £>=2.0. 



L L 


Figure 6. (Color online) The string order parameters = 

(cr|,_icrf, •••cr|j_icr|j) and (L) = (cr|iCr|i+i---crfjCr|j + i) with: (a) J//=l, 

J/=1.2, and D=2 (a point in the odd-Haldane phase); (b) Jh=^, J/=1.6, and D=2 
(a point in the even-Haldane phase). L = \j — i\ is the length of the string, and a 
operator denotes a spin-1/2 Pauli matrix. 
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order parameter ) in the odd-Haldane (even-Haldane) phase decay quickly 

to zero within the lattice distance L=50 (see Fig. ; On the other hand, (^Qz,even 

) in the odd-Haldane (even-Haldane) phase oscillates, whose absolute values converge to 
nonzero values in the large L limit (Fig. [^. Collecting the converged (absolute) values 
of and and show them in Fig. one can clearly see that 01’°'^'^ and 

serve as well-defined order parameters of odd- and even-Haldane phases, respectively. 
In addition, in the vicinity of the critical point, we hnd the (converged) string orders 
obey a power-law scaling versus parameter 5 = \ Ji — Jj\\ 

01 ~ 6 ^“, ( 6 ) 

where (3 = 1/12 [see Fig. (a) and (b)]. 

Subsequently, we analyze the entanglement spectrum (ES) [33], which is defined as 
the -Log 2 (p), where p is the reduced density matrix of half-infinite chain. In the MPS 
framework, the ES can be evaluated through -2 Log 2 A“ (-2 Log 2 A^), corresponding to 
the ES by cutting an odd or even bond. The A“ and A“ denote diagonal matrices on 
odd and even bonds. The ES results of odd (even) bonds are shown in Fig. (a) and 
(b), respectively. Different spectrum structures in the odd- and even-Haldane phases are 
observed. In the odd-Haldane phase, doubly degenerate ES exists on odd bonds (see Fig. 
[^(a)). However, in the even-Haldane phase, the ES becomes doubly degenerated on even 
bonds, instead. It means that the odd- and even-Haldane phases can be characterized by 
doubly degenerate ES on odd and even bonds, respectively. The degeneracy of ES is a 
reminiscence of that in spin-1/2 one-dimensional quantum compass models, where there 
exists a disordered phases with nonzero string orders [Ml [32] . The double degeneracy 
of the ES reflects that the geometric bond space supports a projective representation of 
global symmetries, including the bond centered inversion symmetry na, the rotational 
SO(3) symmetry [35], etc. 


3.3. BLOCK ENTANGLEMENT ENTROPY AND CENTRAL CHARGE 


From the phase diagram shown in Fig. one can find that there exist four 
different phases, and all the QPTs on the phase boundaries separating them are found 
to be critical. In order to study their criticalities, in Figs. ms we show the block 
entanglement entropy Sl, which measures the entanglement between a block of L spins 
and the remaining (infinite) environment. The block entanglements of gapped phases 
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Figure 7. (Color online) Absolute values of and (at L=200), 

which serve as order parameters characterizing the odd- and even-Haldane phases. 



1/6 „l/6 


Figure 8. (Color online) String order parameters (a) amj (b) Q^ieven) ^ 

function of S = Jj — Jf. 
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Figure 9. (Color online) Entanglement spectra on (a) odd bond (-2 Log 2 A“) with 
-D=2.0 and (b) on even bond (-2 Log 2 A^) along the line D=2.0. 



Figure 10. (Color online) Block entanglements (Sl) with Jh = 1 at PI {D=0.7, 
Ji=-2.M), P2(£>=0.2, J/=2.0), and P3 (D=0.536, J/=10) (see PI, P2, and P3 in Fig. 
[^. The block entanglements of these points are fitted well by function k + |Log 2 L, 
and their central charges (c) are determined to be about 1/2. 
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Figure 11. (Color online) Block entanglements (Sl) with Jh = 1 at P4 (19=2.0, 
Ji=lA2) and P5(I9=5.0, J/=1.16). The block entanglement entropies of both points 
are fitted well by fc + |Log 2 L, and their central charges are determined to be c ~ 1. 

are found to be saturated when L is large enough, well satisfying the so-called area 
law [3S]- However, for the critical points, the block entanglement has a logarithmic 
correction to area law, as derived in Ref. 133, Sl increases with L as 

‘S'l =/c + ^^^log2 L, (7) 

where k denotes the intercept and c (c) is the holomorphic (antiholomorphic) central 
charge in the conformal held theory (CFT). The central charges c and c are important 
characterizations of CFT, and can be used to classify the universality class of QPTs 


Five representative points labeled as PI, P2, ■ • P5 (see Fig. are selected from 
the phase boundaries as examples. Among them, PI, P2, and P3 are three points on 
the phase boundaries between two Haldane phases and the AFSP and AF phases, and 
P4, and P5 are critical points between the odd- and even-Haldane phases. The block 


entanglement entropies of PI, P2, and P3 are plotted in Fig. 10, where logarithmic 


divergent behaviors are observed at these three points, and the central charges are 
determined to be c ~ 1/2. It means the QPTs from two Haldane phases to the AFSP 
and AF phases belong to the same universality class, i.e., the Ising universality class. 
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and they can be described by a free fermionic field theory [39], with central charges 
c=c ~ 1/2. On the other hand, logarithmic divergent behaviors are also observed at P 4 


and P 5 (see Fig. 11), but with central charges determined as c ~ 1, showing that the 


topological phase transition between the odd- and even-Haldane phases belongs to the 
Gaussian universality class, and may be described by a free bosonic held theory. 


3.4- Revisit of the HIAC with DM interactions on odd bonds only 

In this subsection, we switch off the DM interaction on even bonds in Hamiltonian 
(|^, and thus reduce the model as 

N/2 

H = . Sh + 

i 

+ ( 8 ) 

As shown in Refs. [221 ESI ES], this model is exactly solvable. Recently Derzhko 
et a/.calculated the ground-state compressibility of a deformable spin-1/2 Heisenberg- 
Ising chain with DM interaction on odd bonds only [lO]. In Fig. 12, we revisit its 


ground-state phase diagram through numerical simulations. Two phase boundaries 
detected by the sharp peaks of bipartite entanglement locate at J/=± 2 a/IP~+”JJ, 
which perfectly agrees with analytical result [21]. As \Ji\ < 2^J+ Jjj, it is a 
disordered phases with nonzero but with vanishing Furthermore, the 

doubly degenerate ES is observed on the odd bonds. These two facts, nonzero 
and doubly degenerate ES, indicate that the disordered phase is an odd-Haldane phase. 
In the region Jj < - 2 /d^T 4 , an AFSP is observed; while an AF phase appears when 
Jj > 2-yiP~+/7j. From the Hamiltonian (js), we realize that the case with positive J/ 
can also be connected to that with negative J/ by a unitary transformation. More 
specihcally, the Hamiltonian with — J/ and that with Jj can be mutually transformed 
by the 7 r-rotation of the spins at sites — 3 and Ai — 2 (or equivalently at — 1 and 
4i) around the x- or y-axis. For instance, if all the spins at sites Ai — 3 and Ai — 2 are 
rotated about x-axis by vr, one can get 

QX _. QtC 

*^ 42 - 3 ( 42 - 2 ) *^ 42 - 3 ( 42 - 2)5 


^4i-3{4i-2) 


S. 


4i-3(4i-2) 


-)■ 


qy 

‘^4i-3(4*-2)> 

- 

‘^4i-3(4*-2)- 


( 9 ) 









J. Phys.: Condens. Matter 


17 



Figure 12. (Color online) Phase diagram of a spin-1/2 HIAC with DM interactions on 
odd bonds [Eq. Q]. It includes three different phases: an AFSP, an AF phase, and an 
odd-Haldane phase. Two second-order critical lines locate exactly at Jf=±2^/D'^ + Jfj 
{Jh is set to be 1 in our calculations). 


After this unitary transformation, the new Hamiltonian is identical to the old one, but 
with — J/. Therefore, when all the spins at sites — 3 and 4i — 2 are rotated over x-axis 

by 71 , the spin conhguration ■■■ — I-1-1- ■ of the AF phase in the region 

J/ > will become an AFSP conhguration ■■•H-|-H-!-■■■ in the 

region Jj < —2^jD‘^ + Jjj, in accordance with our numerical observation. 


3.5. Evolution of the even-Haldane phase 


From the phase diagrams in Fig. and Fig. 12, we hnd that an even-Haldane phase 
will be induced as the DM interactions on even bonds is switched on. In order to explore 
the evolution of such an even-Haldane phase, we introduce the following Hamiltonian 
with alternating DM interactions to interpolate both limiting cases: 

N/2 

H = 'r|J„S2 ._i. Sh + D(Si_,Sl - 






( 10 ) 


The 7 denotes the relevant strength of the DM interactions on even bonds compared 
with that on odd bonds, and thus the values of interest are restricted within the range 
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[0,1]. When 7 = 0 , it reduces to the Hamiltonian (|^ with DM interactions on odd bonds 
only and turns into Hamiltonian (|^ when 7 = 1 . In addition to the studied cases of 7=0 
and 1, we also select 7 = 0 . 2 , 0.3, 0.5 and 0.8 to demonstrate the evolution of the even- 
Haldane phase. The corresponding phase diagrams by tuning D and J/ are provided in 


Figs. 13 (a-d), respectively. Since the effect of the even-bond DM interactions on the 
AFSP is negligible (see Fig.and Fig. 12), we just focus on the region Jj >0. We hnd 
that, an even-Haldane phase can be induced as long as the DM interactions on even 
bonds are taken into account. When 7 is small, the even-Haldane phase locates in the 
region with very large D and Jj (hgures are not shown here). When 7 increases, the 
even-Haldane phase expands quickly while the AF and the odd-Haldane phases shrink. 
As 7 is large enough, this even-Haldane phase comes into the region with small D and 


Jj [see Fig. 13 (a-d)]. Because the even-Haldane phase, the odd-Haldane phase and AF 
phase are adjacent to each other, a tricritical point exists between them. As 7 increases, 
such a tricritical point moves towards the region with small D and J/, and appears in the 


phase diagram with small D and Jj consequently (see Fig. 13 (c) and (d)). It is worth 
noting that the exact position of such a tricritical point is difficult to be determined 
precisely. 


4. DISCUSSIONS AND CONCLUSION 

The ground-state phase diagram and the QPTs in the spin- 1/2 HIAC with uniform 
DM interactions have been investigated by the MPS method. By calculating the odd- 
and even-string order parameters, two kinds of Haldane phases, i.e., odd- and even- 
Haldane phases, have been identified. Furthermore, doubly degenerate entanglement 
spectra on odd and even bonds are observed in odd- and even-Haldane phases, 
respectively. A rich phase diagram including four different phases, i.e., an AFSP, an 
AF phase, odd- and even-Haldane phases, have been obtained. The TQPT between 
the odd- and even-Haldane phases is with central charge c=l (a Gaussian type phase 
transition) and with critical exponent (3 = 1/12. The central charges on the other phase 
boundaries are c=l/2, therefore the QPTs from nonmagnetic (two Haldane) phases to 
the magnetic ordered (AFSP and AF) phases belong to the Ising universality class. 

In addition, the bipartite entanglement entropy has been shown as a very powerful 
tool for capturing QPTs (including TQPTs). Two kinds of nonlocal string orders and the 
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Figure 13. (Color online) Phase diagrams of a spin-1/2 HIAC with DM interactions 
on odd and even bonds for (a) 7 = 0 . 2 ; (b) 7 = 0 . 3 ; (c) 7 = 0.5 and (d) 7 = 0 . 8 . The 7 
denotes the relevant strength of the DM interactions on even bonds compared with 
that on odd bonds [see definition in Eq. @]. 

doubly degenerate ES can be used to distinguish the odd- and even-Haldane phases. The 
nonzero string order parameters imply the breaking of the hidden topological symmetry 
in the odd- and even-Haldane phases. It is worth noticing that the existence of the 
odd-Haldane phase has been already proved for the Heisenberg-Ising chain without 
Dzyaloshinskii-Moriya interaction H. We would like to remark that the tricritical 
point (see the hlled square ■ in Fig. is rather difficult to perform an accurate 
calculation, posing a challenging problem in numerical simulations. In order to explore 
the evolution of the even-Haldane phase, we introduced a Hamiltonian with a tunable 
DM interactions on even bonds. Our results show that the even-Haldane phase can be 
induced so long as the DM interactions on even bonds are taken into account. The DM 
interactions on even bonds plays a key role in the formation of the even-Haldane phase. 

Above we only considered the case with positive DM interaction (D > 0). In fact, 
we realize that the case with negative D is connected to that with positive D via a 
transformation: one can rotate all the spins by vr angles about the x-axis, he., —)■ 

—)■ —S^ and —>• —S^. This rotation leaves the Hamiltonian intact, but adds a 

minus sign before the DM term. Therefore, it is sufficient to consider only the case of 
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